Charge transfer, confinement, and ferromagnetism in LaMnOs/LaNiOs 

(OOl)-superlattice 
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Using first-principles density functional theory calculations, we investigated the electronic struc- 
ture and magnetic properties of (LaMn03)m/(LaNi03)n superlattices stacked along (OOl)-direction. 
The electrons are transferred from Mn to Ni, and the magnetic moments are induced at Ni sites that 
are paramagnetic in bulk and other types of superlattices. The size of induced moment is linearly 
proportional to the amount of transferred electrons, but it is larger than the net charge transfer 
because the spin and orbital directions play important roles and complicate the transfer process. 
The charge transfer and magnetic properties of the {m^n) superlattice can be controlled by changing 
the m/n ratio. Considering the ferromagnetic couplings between Mn and Ni spins and the charge 
transfer characteristic, we propose the (2,1) superlattice as the largest moment superlattice carrying 
^8/jb per formula unit. 

PACS numbers: 75.70.Cn, 73.20.-r, 75.47.Lx, 71.15.Mb 



Introduction — Recent advances in the layer-by-layer 
growth technique of transition metal oxide (TMO) het- 
erostructures have created considerable research interest 
[l|,[2|. In TMO, multiple degrees of freedom (i.e., charge, 
spin, orbital) are coupled to each other, often creating 
novel material characteristics such as high-temperature 
superconductivity and colossal magneto resistance [3|- 
By making artificial heterostructures of TMO, one can 
control those degrees of freedom and band structures, 
and therefore create or design the new 'correlated elec- 
tron' properties. Previous TMO superlattice studies [j- 
0] have shown that various unexpected material phenom- 
ena can be realized at the TMO heterointerface such as 
magnetism |9l-[ll| and superconductivity p^ . 



In this context, the superlattices composed of LaNiOs 
(LNO) and LaMnOs (LMO) are of particular interest. 
A recent experiment by Gibert and co-workers reported 
the exchange bias in LMO/LNO stacked along the (111)- 
direction [13|- Mn-to-Ni charge transfer is expected at 
the LMO/LNO interface, which may cause a sizable mag- 
netic moment in the Ni ions even if LNO is paramagnetic 
in bulk [3| and many other superlattices |1#J16| . Impor- 
tantly, however, it is quite unclear if the same mechanism 
is also working in the LMO/LNO superlattice stacked 
along (001). Although the main concern of the paper 
by Gibert et al. is the (lll)-structure, their data does 
not seem to support the same physics taking place in the 
(OOl)-stacked LMO/LNO. On the other hand, a recent 
extensive experimental study by Hoffman et al. [I3 re- 
ports that the same type of charge transfer also occurs 
in the (OOl)-case, and the magnetic signals were clearly 
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observed from Ni sites. Furthermore, a recent theoret- 
ical work by Dong and Dagotto [18| suggests a differ- 
ent mechanism for the induced magnetic moment (M) 
at Ni. Their tight-binding calculations show that the in- 
duced magnetic moment in the (lll)-superlattice is bet- 
ter understood by the quantum confinement effect rather 
than by the charge transfer. While the effect of confine- 
ment is strongest in the (111)- and weakest in the (001)- 
structure, this study also raises an important question re- 
garding the induced Ni moment in the (OOl)-superlattice. 
However, due to the lack of first-principles calculations 
for the (OOl)-structure, the detailed understanding of this 
system has not yet been achieved. 

In this paper, we examine the (OOl)-superlattice with 
first-principles density functional theory calculations. 
Our calculations of (LMO)m/(LNO)n with several com- 
binations of (m^n) clearly show that the significant charge 
transfer occurs and the magnetic moments are induced at 
Ni sites as in the (lll)-case. Furthermore, the size of an 
induced moment is linearly proportional to the amount 
of charge transfer. However, the simple count of net elec- 
tron transfers cannot explain the size of the moment be- 
cause the transfer process occurs in a complicated way 
that depends on spin and orbital directions. The ma- 
jority spin and d^z'^_j.2 orbital are the major channels in 
the electron transfer while the occupation changes in the 
minority spin and 6^3.2 _^2 orbital are not negligible. We 
also found strong ferromagnetic (FM) couplings between 
Ni and Mn, whereas the Mn-Mn and Ni-Ni spins are an- 
tiferromagnetically aligned in some cases. As a result, 
superlattices with (77i,n) = (l,l), (1,2), and (2,1) are al- 
ways FM regardless of U value, which can have important 
implications for application. Considering the amount of 
charge transfers and FM couplings across the interface, 
the (2,1) structure is proposed to have a largest moment 



of ^ 8/J.B' 

Computation Details — For calculating 

(LMO)m(LNO)n (OOl)-superlattices (2 ^ m ^ 3, 
2 ^ n ^ 4), we used the projector augmented wave 
(PAW) potentials [19| and generalized gradient ap- 
proximation (GGA) proposed by Perdew [20| for the 
exchange-correlation functional, as implemented in 
the VASP code [21]. To study the effect of electron 
correlation, we also used the GGA-\-U scheme within 
the rotationally invariant formalism and the double- 
counting formula, as firstly proposed by Liechtenstein 
et al. [22|. We used four different sets of U and J 
for La-4/, Ni-3(i, and Mn-3d states: i) /7aii = Jaii 
= ii) /7La=6eV,JLa=0.5eV, t/Ni=6eV, JNi=0.5eV, 
^Mn=5eV, JMn=0.5eV, iii) /7La=3eV, JLa=0.5eV, 
^Ni=3eV, JNi=0.5eV, /7Mn=2.5eV, JMn=0.5eV, and iv) 
f/La=OeV, JLa=OeV, /7Ni=leV, JNi=OeV, /7Mn=4eV, 
'>^Mn=leV [23|. Note that the last setup for U and J 
is the one used by Gibert et al. for the (lll)-structure 
|13j . While we are mainly presenting and discussing 
the results from U^i=0 and 3 eV, it was found that 
the main claims and conclusions are not changed in 
the other sets of parameters. The wave functions were 
expanded in plane waves with a kinetic energy cutoff of 
500 eV. We used a k-point set generated by the 8x8x4 
Monkhorst-Pack mesh for the (1,1) superlattice and 
used equivalent k-points for other (m, n)-superlattices. 
Ionic coordinates were optimized until the residual forces 
were less than 0.01 eV/A. Wigner-Seitz radii of 1.286 
and 1.323 A were used for the projection of Ni and Mn 
atoms, respectively, as implemented in the VASP-PAW 
pseudopotential. We assumed that the LMO/LNO 
superlattice is grown on the SrTiOa substrate by setting 
the in-plane lattice constant fixed at a=6=3.905A. We 
used the tetragonal supercell and the optimized c-lattice 
parameters for each (m^n) superlattice within the FM 
spin configuration. 

Bulk and structural property — The bulk LNO is known 
to have low-spin d^ electronic configuration and to remain 
as a paramagnetic (PM) metal down to low temperature 
Q. The local density approximation (LDA) and GGA 
calculation (/7=0) predict the correct PM ground state 
for the bulk phase and some other superlattice struc- 
tures such as LNO/LAO and LNO/STO Blllj, while 
LDA+t/ predicts the local moment formation at Ni site 
[H, [23. In our calculations, GGA+/7 yields the Ni mo- 
ment of 1.10 and 1.36 fiB for U=3 and 6eV, respectively. 
In bulk LMO, Mn^+ has high-spin d"^ configuration in 
which is, tlgO^^. In GGA {U=0) calculation, it is found 

that the small amount of ej electron is transferred to t2g 
state due to the down-spin t2g bands close to the Fermi 
level. The calculated magnetic moment is increased from 
3.46 /iB at U=0 to 4.05/iB at U=6 eV. 

The optimized out-of-plane lattice parameter of bulk 
LNO and bulk LMO (with the fixed a, 6 lattice of STO 
value) are found to be clno = 3. 798 A and clmo = 
3.918A, respectively. In the (LMO)m/(LNO)n superlat- 
tice, the Ni-O-Ni distance (<iNi-Ni) and the Mn-O-Mn dis- 
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FIG. 1: (Color online) (a) Atomic structure of 

(LMO) 2 /(LNO) 2 superlattice. Grey, purple, green, and red 
colors stand for Ni, Mn, La, and O atoms, respectively. 
c/tmi-tm2 denotes the distance between the TMl-plane and 
TM2-plane, where TMl and TM2 are Mn or Ni atoms, (b) 
The calculated ^Mn-Mn, c/Mn-Ni, and c/ni-ni distances. The 
values for bulk LNO and LMO are indicated by (??i,n) = (0,l) 
and (1,0), respectively. Note that in (2,4), two different types 
of c^Ni-Ni exist. The shorter g^ni-ni corresponds to the Ni-Ni 
distance between the two inner-most layers. 



tance (<iMn-Mnin Fig- 1(a)) along the c-axis are changed 
so that Clno < c^Ni-Ni and clmo > <^Mn-Mn- As a result, 
the distances between the two transition metals (TMs) in 
the superlattice are always larger than clno and smaller 
than Clmo, as clearly shown in Fig. [TJb) for the U=0 
case. It is noted that the inner layer (iNi-Ni approaches 
to Clno as the thickness of LNO layers, n, increases. We 
also found the same trend in the U > results. 

It is found that TM-O-TM bond angles in 
(LMO)^/(LNO)n are generally not 180° [H. The in- 
plane angle between Mn-O-Mn (ZMn-O-Mn) in (?72=l,n) 
superlattices is ^180° since these superlattices have mir- 
ror symmetry with respect to the Mn02-plane. On the 
other hand, for (771=2, n) structure, ZMn-O-Mn decreases 
as n increases. Similarly, the in-plane angle between Ni- 
O-Ni (ZNi-O-Ni) at the interface of (m=2,n) superlat- 
tices also decreases as n increases, while ZNi-O-Ni?^ 180° 
in (l,n) superlattices. It is noted that ZNi-O-Ni is in- 
creased for the bulk-like Ni atoms. To see the change 
of the out-of-plane TM-O-TM bond angle and the pos- 
sible octahedra rotations, we performed the geometrical 
optimizations from distorted structures as starting ge- 
ometries in which the atomic positions are shifted to- 
ward the in-plane oxygens (with no change in the lattice 
parameters). It was found that the O atoms return to 
their original position and the out-of-plane bond angles 
between Ni-O-Ni remain as 180°. 

Charge transfer and Ni magnetic moment — Figs.[2fa)- 
(b) summarize the calculated result of the charge trans- 
fer between Ni and Mn for several (m^n) combinations 
of (LMO)m/(LNO)n, where up and down panels corre- 
spond to the gain and loss of electrons, respectively. The 
number of TM-d electrons in the bulk LNO and the bulk 
LMO are set to be zero as a reference point for Ni and 
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FIG. 2: (Color online) The amount of charge transfer in 
(LMO)m/(LNO)n superlattices (denoted by {m,n)) calcu- 
lated with (a) [7=0 and (b) U—3eV. Red boxes and blue 
circles represent the calculated charges of Mn and Ni atoms, 
respectively. Zero-charge indicates the values from the bulk 
LMO and LNO. (c)-(d) The amount of transferred charge to 
Ni as a function of U. (2,3) s and (2,4) b stand for the bulk- 
like Ni atom in the (2,3) and (2,4) superlattices, respectively, 
while others (with no subscript) refer to the interface Ni. 



Mn charge, respectively. The results correspond to the 
most stable spin configuration among all possible spin or- 
ders for given (m^n) structures. A clear common feature 
is that the electrons are transferred from Mn to Ni. Al- 
though the amount of charge transfer in the (lll)-case is 
not given in Ref. [l3|, we expect that the charge transfer 
in the (lll)-superlattice is larger than that of the (001)- 
case because the (lll)-interface creates more Mn-O-Ni 
bonds than the (OOl)-interface does. This point is also re- 
flected in the result of the magnetic moment, which will 
be discussed further. The transferred electrons mostly 
reside at the interface Ni sites and the valence change in 
the bulk-like (inner layer) Ni is relatively small, as clearly 
seen in the result of (2,3) and (2,4). 

Note that, since Mn donates electrons to Ni, the 
amount of charge transfer and the Ni valency can be con- 
trolled by changing the superlattice composition {m^n). 
For a larger ratio of m/n^ the induced change in the Ni 
valence becomes larger while that for the smaller m/n the 
change smaller. By comparing (1,1) structure with (2,1), 
one can flnd that the Ni-d occupation is larger in (2,1), 
where the two Mn ions can provide electrons to one Ni. 
The same feature is conflrmed by comparing (2,2) with 
(2,3). As we will discuss below, the transferred electrons 
induce the magnetic moment at Ni, and therefore, the 
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FIG. 3: (Color online) (a)-(b) The calculated magnetic mo- 
ments of Mn and Ni in the superlattices with (a) U—^ and (b) 
[/=3eV. Red boxes and blue circles represent the Mn and Ni 
values, respectively, (c)-(e) The calculated Ni magnetic mo- 
ment as a function of transferred charge to Ni for (c) [/=OeV, 
(d) [/=3eV, and (e) ?7=6eV. Filled symbols represent the in- 
terfacial atoms, whereas open symbols represent the bulk-like 
ones 



magnetism can also be controlled by changing the super- 
lattice compositions {m^n). 

The main feature regarding the charge transfer is main- 
tained even when the on-site correlation U is turned on 
as shown in Fig. [2fb). The same curve shapes are found 
as in the U=0 results (Fig. Efa)), indicating the same 
type of charge transfer. The effect of U is to reduce the 
amount of charge transfer. /7=6 eV results [23] are also 
found to be consistent with U=3 (Fig. [2fb)). The effect 
of correlations that reduces the charge transfer is more 
clearly seen in Figs. [2fc) and (d), where the increase of 
Ni-d occupation (with respect to the bulk value) is plot- 
ted as a function of U. The decreasing feature is evident 
for all compositions of (m^n) and for both interfacial and 
bulk-like Ni. It is noted that, in the bulk-like Ni sites, 
the valence change caused by charge transfer is close to 
zero if U is large enough (see the dashed line in Fig. [2fd) 
at [/=6 eV). 

For the (lll)-superlattices of LMO/LNO [13], it is re- 
ported that the magnetic moment is induced at Ni, which 
is originally paramagnetic in bulk, and the exchange bias 
is manifested by this induced moment. For the (001)- 
case, however, it is unclear if the Ni atoms are spin po- 
larized. It seems that two experimental studies arrive at 
the different conclusions regarding this question for the 
(OOl)-case. Gibert et al. [13| reported that the exchange 
bias is not manifested in the case of the (OOl)-interface, 
while Hoffman et al. [I3 clearly observed the magnetic 
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FIG. 4: (Color online) The calculated charge transfer for the 
minority- spin {U=Oe'V). Red boxes and blue circles represent 
the Mn and Ni values, respectively. Note that the charge 
transfer shape is opposite of that of Fig. 2(a). 



signals. Therefore, a detailed theoretical analysis is re- 
quired to understand the magnetic property of the (001)- 
structures. 

Our calculations clearly show that the Ni magnetic mo- 
ment is also induced in (OOl)-superlattices. After calcu- 
lating all the possible spin orders for given (7Ti,n), we 
present the most stable spin configurations in Figs. [3fa)- 
(b), where up and down panels represent majority (up) 
and minority (down) spin, respectively. It is noted that 
the Ni ions have non-zero spin moment even in the U=0 
calculations (Fig. [3ja)). 

The calculated Ni moment is ~0. 08-0. 51 fiB at [/=0, 
and enhanced up to ~ 1.10-1. 53 /j^b at U=3 eV. It is noted 
that, for the (2,2) superlattice, the calculated value of 
Mni is 0.47 jiiB at [/=0, similar to the experiment ~0.35 
A^B [i3- ^^ fo^ ^Mn , there is a significant difference 
between the calculated value of 3.14/iB and the exper- 
imental one {^2iiB ll3)- Although the origin of this 
discrepancy is unclear, we emphasize that the calculated 
value is in good agreement with Mn charge status of 4+, 
which is supported both by our calculation (see Fig.^ 
and the x-ray absorption spectroscopy data in Ref. [l3 
|23j . It should be noted that for bulk LNO and other 
superlattices such as LNO/LAO 0, IH and LNO/STO 
[iq . GGA (or LDA; U=Q) predicts zero moment for Ni. 
Therefore, our result of finite Ni moments at /7=0 is a 
clear evidence for the induced net moment. 

It is instructive to compare the magnetic property of 
the (OOl)-superlattice with (111). In both cases, Ni mag- 
netic moment is induced and coupled to Mn spins ferro- 
magnetically. Also, Mn-Mn and Ni-Ni coupling is anti- 
ferromagnetic (AFM) in (2,2) structure. The notable dif- 
ferences are found in the size of the magnetic moments. 
According to Gibert et al. [H, 1.1 < Mni < IAjib for 
the (lll)-case. These values are much larger than our 
results for the (OOl)-interface with [/=0 (see Fig. [3fa)). 
To be more precise, we performed the calculations with 
the same U and J values as used in Ref. [l3|. The result 
clearly shows that the calculated Ni moment is always 
smaller in (OOl)-superlattice than in (111) by ^ O.Sjub^ 
while the Mn moment is larger in (001) by ~ 0.3/iB- Note 
that this trend is also consistent with the charge transfer 
feature, as discussed above. 



The origin of the induced moment in LMO/LNO is un- 
der debate [isl, [l3, [l8| . Gibert et al. speculated about 
the charge transfer and two-dimensional confinement as 
a possible origin of the induced Ni moment and seem 
to have concluded that neither of them plays a signifi- 
cant role [13|- On the other hand, in an interesting re- 
cent study, Dong and Dagotto suggest that the induced 
moment is better understood as a result of the spin- 
dependent quantum confinement rather than the charge 
transfer especially for the case of the (lll)-superlattice. 
This confinement effect is shown to be strongest in the 
(lll)-interface and weakest in (001). Thus, further study 
seems necessary for the (lll)-case, and it is important to 
understand the role of confinement and charge transfer 
in the (OOl)-case. 

To address this point, we present the induced Ni mo- 
ment as a function of the amount of electron transfer in 
Figs. [3jc)-(e). A linear dependency is quite clear, es- 
pecially for the non-zero U calculations, and the U=0 
result is not very far from the linear fit. This point sug- 
gests that the charge transfer is the key origin of the 
induced Ni moments in case of the (OOl)-superlattices. 
It is noted that while the size of Ni moments is larger for 
/7 > 0, the moment enhancement by heterostructuring 
is larger in U=0 because the correlation U reduces the 
charge transfer. It is also consistent with the picture of 
charge-transfer-driven moment formation. 

Spin and orbital dependency — It is interesting to note 
that the charge transfer between Mn and Ni takes place 
in a different way in the minority-spin bands from the 
majority. It is found that the number of minority-spin 
electrons at Mn is enhanced in the superlattice (com- 
pared to the bulk value) whereas that of Ni is mostly 
reduced, which is opposite to the case of majority-spin. 
Fig. m^a) clearly shows that the sign of charge transfer is 
reversed in the minority spin case (compare Fig. |4] with 
Fig. [Sfa)). As a result, the induced Ni moment is larger 
than the net charge transfer, and we note that this ef- 
fect is largest for the (1,1) superlattice with the largest 
induced Ni moment. 

Importantly, the total amount of charge transfer is 
dominated by majority-spin channels while the change 
in the minority-spin electron occupations is relatively 
small. We note that this point is consistent with the 
spin-dependent quantum confinement picture suggested 
by Dong and Dagotto [18|, even though their analysis 
is best applied to the (lll)-interface, and the effect is 
relatively weak in the (OOl)-case. In this picture, the 
majority-spin Ni-e^ wavefunction is more widely spread 
out while the minority-spin electron is localized. The 
delocalized feature of the majority-spin bands and the 
more overlap with the neighboring up-spin Mn bands are 
consistent with our results that the charge transfer oc- 
curs through the majority-spin channel. For the major- 
ity spins, the Ni-d occupations decrease as U increases. 
In the minority spin bands, the occupation changes are 
much smaller [23| . Although the occupation changes are 
quite small, interestingly the amount of occupation en- 



hancement is found to increase as U increases, which is 
an opposite trend to the majority spin case. 

By analysing the orbital occupations, we found that 
the major charge transfer channel is d^z^-r^ orbital. 
Compared to the bulk value, the majority-spin Mn- 
d^z'^_r2 occupation is reduced by ~0. 07-0. 16 depending 
on (m^n) {U=0 eV). The change in Mn-dj;2_y2 occupa- 
tions is less reduced as ~ 0.02-0.11. Due to its wave- 
function shape, d^z'^_r2 can have more overlaps and be a 
more efficient channel for this process. The same feature 
is found in the change of Ni-e^ occupations. Ni-d^z'^-r'^ 
occupations get enhanced by ~ 0.12-0.18 while the Ni- 
dx2_y2 by ~ 0.06-0.10. The main features of this spin and 
orbital dependent charge transfer are maintained also in 
the GGA^U calculations [23^. 

Designing ferromagnetic superlattices — Making an 
FM TMO superlattice with a large magnetic moment 
and a high Curie temperature is an important issue for 
applications [25n28]. It is noted in our system that the 
induced Ni moment can ferromagnetically align with Mn 
spins as in the (lll)-superlattice of LMO/LNO [H. As 
summarized in Figs.[3fa)-(b), our calculations show that 
in the (OOl)-superlattice FM coupling across the inter- 
face {i.e., between Ni and Mn) is always favored ener- 
getically. That is, the interface FM spin arrangements 
always have less total energy than AFM ones, regardless 
of the other parts of spin orders and independent of U 
values (see Figs. [3fa)-(b)). For example, the FM (1,1) 
superlattice has the lower total energy than the AFM 
one by 130 meV/(LNO)i(LMO)i at /7=3eV, which cor- 
responds to the magnetic coupling JNi-Mn= 58 meV with 
5'Mn=l-6 and 5'Ni=0.7. While the Ni-Mn spins are al- 
ways aligned ferromagnetically, the Mn-Mn and Ni-Ni 
couplings are either AFM or FM depending on U and 
{m,n) (see Figs. [3i;a)-(b)) [231]. 

Our result has an interesting implication in regard to 



the design of the magnetism of superlattices. Since the 
interface (Ni-Mn) coupling is always FM, the superlat- 
tice compositions of (1,1), (1,2), and (2,1) should be FM 
carrying large total moments. Furthermore, (2,1) struc- 
ture is expected to have the largest moment, especially 
because the amount of charge transfer will be largest in 
this case as we discussed already (that is, largest m/n ra- 
tio). Therefore, based on our analysis for charge transfer 
and magnetic coupling, one can expect the (2,1) structure 
to be the large moment FM superlattice. The calculated 
total moments of (1,1), (2,1), and (1,2) are 4.14, 7.71, 
and 4.14/iB, respectively, at U=0 eV. It is noted that 
the significant amount of magnetic moment can actually 
be controlled by changing the (m^n) compositions. The 
calculated moments by U=3 are 5.0, 9.0, and 6.0/iB for 
(1,1), (2,1) and (1,2), respectively, being consistent with 
U=0 results. The suggested guiding principle based on 
the FM coupling and the charge transfer can be useful to 
design the magnetic superlattices. 

Summary — Our ffist-principles calculations show that 
the magnetic moments are induced at Ni atoms in the 
(OOl)-oriented (LMO)m/(LNO)n. The induced Ni mo- 
ment is governed by the electron transfer from Mn to Ni 
and the amount of charge transfer increases as m/n in- 
creases. Spin and orbital directions also play important 
role. Our analysis based on the FM couplings between 
Mn and Ni and the charge transfer features can provide 
an useful designing principle for the magnetic TMO su- 
perlattices. 
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